Large scale numerical simulations of "ultrametric" long-range depinning 
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The depinning of an elastic line interacting with a quenched disorder is studied for long range 
interactions, applicable to crack propagation or wetting. An ultra-metric distance is introduced 
instead of the Euclidean distance, allowing for a drastic reduction of the numerical complexity 
of the problem. Based on large scale simulations, two to three orders of magnitude larger than 
previously considered, we obtain a very precise determination of critical exponents which are shown 
to be indistinguishable from their Euclidean metric counterparts. Moreover the scaling functions are 
shown to be unchanged. The choice of an ultrametric distance thus does not affect the universality 
class of the depinning transition and opens the way to an analytic real space renormalization group 
approach. 



I. INTRODUCTION 

The depinning of an elastic interface in a random en- 
vironment gives a common theoretical framework to de- 
scribe physical phenomena as various as the advance of a 
magnetic wall, the propagation of a fracture front or the 
wetting of a disordered surface (see e.g. Ref. jH 0- Q 
for a recent review). The richness of the physics en- 
countered in these different situations results from the 
same feature. The disorder of the environment which 
tends to anchor the interface competes with an elastic- 
like term which tends to keep the interface smooth. The 
tuning of an external driving force allows to go through 
a critical transition. Below threshold the interface can 
only advance over a finite distance before stopping in a 
blocked conformation. Above threshold the interface can 
move freely and acquire a finite velocity. At threshold 
the system is characterized by a set of universal critical 
exponents. 

In the case of over-damped dynamics the motion can 
be described by the following stochastic equation: 

lidth(x,t) = F ext (t) + F e i(x,t) +j[x,h(x,t)] (1) 

where h(x,t) denotes the position of the front, F ext (t) 
the external driving force, F e i(x,t) the elastic term due 
to the distortion of the front and j[x,h(x,t)] the frozen 
disorder. Depending on the physical phenomenon con- 
sidered the elastic term can take different forms. In the 
case of magnetic walls or of wetting in a Hele-Shaw cell, 
the interactions to be taken into account are short ranged 
and at first order the elastic force can be estimated by a 
simple Laplacian term: F e i(x, t) = V 2 h(x, t). In the case 
of the advance or receding of a triple line on a disordered 
substrate Q or of the propagation of fracture front long 
range interactions have to be considered. In the latter 
case elastic interactions are mediated via the bulk along 
the whole front. A first order perturbation analysis for 
the front roughness gives 
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This long range elastic string model has been widely 
studied over the last ten years. In particular the interface 
is shown to exhibit a self-affine roughness: the width w 
of the interface scales with the system size L as w oc 
The different numerical works |fj|, IjJ, |a |9( performed give 
estimations of the roughness exponent ( spreading in the 
interval [0.34 — 0.40]. The latter exponent and more gen- 
erally the universality class of the model depend strongly 
on the long range nature of the kernel in Eq. (2). For 
instance, if F e i oc V 2 h, then £ = 1.25. The most recent 
results obtained by Rosso and Krauth Q give a value 
C ~ 0.39 significantly larger than the theoretical predic- 
tion £ = 1/3 obtained by one- loop calculations of renor- 
malization group technique [Tol 111) and equally smaller 
than the recent two-loop estimation^] C ~ 0.47. Note 
however that these values are not consistent with the re- 
sults ( ~ 0.5 — 0.6 obtained experimentally for interfacial 
fracture0 or wetting [l3l H^ . Therefore it is of utter- 
most importance to have an accurate determination of 
these critical exponents, and thus to be able to study 
large system sizes. 

In the following we present an ultrametric version of an 
extremal model of depinning (see Ref. y| for a detailed 
study of the original model with Euclidean metric) . The 
complexity of the elementary step of calculation is shown 
to scale with the system size L as log 2 L instead of L in 
the original model. This allows us to perform simulations 
on systems of size L — 2 20 rj 10 6 , which corresponds to 
a gain of two to three orders of magnitude compared 
with other published works based on the Euclidean met- 
rics. The universality class of the model is shown to be 
unchanged. Beyond the numerical acceleration, this ul- 
trametric model of depinning may thus also serve as a 
starting point for a real space renormalization analysis. 

The paper is organized as follows: In a first part we 
recall the definition of the original model, we then define 
the ultrametric version and present the main features of 
the new model. In the second part we give results of 
simulations performed on large systems and focus on the 
critical properties of this ultrametric model of depinning. 
We conclude that it lies in the same universality class 
as its euclidean version. The details of the numerical 
implementation of the algorithm are finally given in an 
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appendix. 



II. EXTREMAL DRIVING OF A DEPINNING 
FRONT 

Various numerical techniques can be used to study 
the depinning phenomenon in the vicinity of the critical 
threshold. Early works used a direct integration of the 
equation above waEulerQ or Runge-Kutta [l5l | schemes. 
Recently Rosso and Krauth developed an iterative al- 
gorithm to determine the blocked conformations corre- 
sponding to a given constant forcing. They used peri- 
odic boundary conditions both along the front and in 
the direction of propagation and the critical threshold is 
reached when the last blocked conformation starts mov- 
ing. 

In the case of strong pinning the advance of the front 
proceeds by successive local instabilities. This avalanche 
behavior is characteristic of the motion of a depinning 
front. This property can be exploited to develop an effi- 
cient algorithm describing the motion of the front close 
to the critical threshold. Instead of driving the front at a 
constant external force, it consists of tuning the latter at 
the exact value such that one and only one site can depin 
at a time. So doing, the sequence of depinning events is 
preserved. The implementation of this extremal dynam- 
ics which has been used since 1992 in various interface 
growth models [H El El E E3, ^ is straightforward. 
At each iteration step one needs to identify the weakest 
site, to advance it up to the next trap and to update 
the long range elastic forces due to the change of front 
conformation. The latter operation scales with the size 
of the front. The great advantage of this method is that 
the system remain constantly at the edge of the critical 
behavior and it is not necessary to tune the external driv- 
ing force (see e.g. Ref. [22.!] for a discussion on the use of 
extremal dynamics to reach the critical state and more 
generally on the link between Self Organized Criticality 
and classical critical transitions). 

Based on the above described extremal modeling, we 
now detail the way to implement the model numerically. 
The front is discretized along a regular horizontal grid; 
i £ [0, L — 1] and hi are the coordinates along the front 
and in the direction of propagation respectively. Traps 
of random depths ji are randomly distributed along the 
direction of propagation. The distortion of the front in- 
duces elastic forces ff via a Green function Gij = G(rij) 
where ry is the euclidean distance separating two sites i 
and j along the front. In the case of a fracture front, a 
discretized version of the elastic redistribution function 
is such that 
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Let us consider a given conformation of the front. For 
each site (i 7 hi) located in a trap of depth ji, we can de- 
fine a local depinning threshold Si = 7, — /• : this site 



depins as soon as the external driving force F overcomes 
the threshold F > Si. The depinning threshold s(t) of 
the front conformation obtained at iteration t thus corre- 
sponds to the minimal external force to be exerted so that 
at least one site of the front can depin: s(t) — min^ Si. Fi- 
nally the critical threshold s* above which the front can 
freely propagate is s* = max t s(t). The basic rule of the 
extremal driving consists simply at each iteration step t 
of simply tuning the external force at exactly the value 
of the depinning threshold of the current front confor- 
mation: F(t) — s(t). Once identified, the extremal site 
is advanced up to the next trap, the elastic forces are 
updated to take into account this local displacement and 
the new value s(t + 1) of the front depinning threshold 
is evaluated. We summarize below the elementary steps 
of the algorithm used to run the model and we estimate 
their complexity respectively to the size L of the system. 



A initialization hi 
i e [0 - L] , L 

B identification of the extremal site io 
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C advance of the extremal site 
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D update of the trap depth ji <- 

E update of the elastic forces ff 
where Gy oc \i — j\~ 2 

F back to step B 



rnd [1] 
-/f + G ii0 6h i0 [L] 



where rnd stands for a random number and a <— b for 
the assignment of value b to variable a. Except the first 
initialization step, the two limiting steps are the identifi- 
cation of the extremal site and the update of elastic forces 
along the front which both scale linearly with the system 
size L. This sequence of elementary steps is then iterated 
T times to obtain statistical averages of the quantities of 
interest. 



III. ULTRAMETRIC DEPINNING 

We now turn to the presentation of the ultrametric 
model. The basic rules of the extremal model remain 
identical but the redistribution of the elastic forces. In- 
stead of using the natural euclidean distance along the 
front we use an ultrametric distance. The structure of 
the algorithm stays roughly similar to the previous one 
but steps B and E are shown to be characterized by a 
complexity in log 2 L instead of L. 

The most natural structure to be used in the context 
of a model with ultrametric distance is a dyadic tree. 
Let us first build such a tree whose final leaves are the L 
sites of the depinning front. As illustrated on Fig. 2] the 
simplest ultrametric distance between two sites i and j is 
the number of branches that composes the shortest path 
on the tree between the two sites. This is exactly twice 
the depth m of the nearest common ancestor of these two 
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FIG. 1: The depth m of the first common ancestor of sites ii 
and 12 is used to compute the ultrametric distance between 
points ii and ii as d u (ii, 82) = 2 m — 1 



sites. Different choices of distance are possible based on 
the tree structure. In the following we use a definition 
that preserves the scaling of the Euclidean distance: 



d u (i,j) = 2 m -l (4) 

An important characteristic of this ultrametric distance 
is its degeneracy. Namely there is one point at distance 
d — 2 1 — 1 = 1, two points at distance d = 2 2 — I = 3 and 
2 P_1 points at distance d — 2 P — 1. For a set of L = 2™ 
points one thus counts only n = log 2 L different values 
for the distance between two points of the set. There lies 
the main advantage of the choice of an ultrametric dis- 
tance from the computational complexity point of view. 
The expression of the elastic Green function then derives 
directly from the definition of the new distance: 

=GK(i,j)]«— i-^ (5) 

<*U \* ! J ) 

log 2 L 

G u = - J2 2 m ~ 1 G [2 m - 1] 

m— 1 

Using this definition, we can easily accelerate the up- 
date of elastic forces (step E). As illustrated on Fig. |2 
instead of updating L sites, we can update only log 2 L 
subtrees corresponding to sites located at the same ul- 
trametric distance of the extremal site. A similar gain 
can be obtained on step B for the determination of the 
extremal site. Technical details regarding the numerical 
implementation of the algorithm are developed in the ap- 
pendix. The basic steps of the algorithm then scale with 
log 2 L instead of L. The price to pay for this numerical 
advantage is the loss of the translational invariance. We 
show below that it does not affect the universality class 
of the model. 

The simulations have been performed on systems of 
sizes up to L = 2 20 . The numerical runs were performed 
over a large duration T respective to the natural corre- 
lation time of the system tl oc L z . In the case of the 
largest system L = 2 20 we used T = 2.510 10 » 100r L . 



FIG. 2: When redistributing elastic forces after the depinning 
of the extremal site io , all sites located at the same ultramet- 
ric distance m = d u (i, io) belong to a common subtree and 
receive the same contribution. This allows an update by block 
that scales with the number log 2 L of these subtrees. Simi- 
larly the tree structure allows to find the next extremal sites 
in only log 2 L operations. 



IV. CHARACTERIZATION OF THE CRITICAL 
STATE: NUMERICAL RESULTS 

The propagation of depinning fronts at the critical 
threshold exhibits a rich phenomenology. The front 
presents a self-affine roughness: its width w(Ax) mea- 
sured over a distance Ax scales as w(Ax) oc Ax 1 ' , where 
£ is the roughness exponent. 

The dynamics of depinning is characterized by an 
avalanche behavior. In the framework of an extremal 
dynamics this can be described by the distribution of 
the distances r between two successive depinning sites: 
P(r) oc r~ a . Generalizing this distribution for sites cor- 
responding to depinning events separated by a given time 
lag allows to obtain in addition the dynamic exponent 
z which characterizes the spreading of the avalanches. 
Namely the lateral extension £ of an avalanche of dura- 
tion At scales as £ oc At 1 ^ . 

Another quantity of interest is the external force s 
needed to depin a given conformation of the front. It 
can be shown that the distribution Q(s) of these 

front depinning forces exhibits a singular behavior close 
to the critical threshold s*: Q(s) oc (s* - s) M . 

In the following we present simulations of ultrametric 
depinning performed on large systems (up to L = 2 20 ). 
We recover all critical features described above with ex- 
ponents numerically indistinguishable from their coun- 
terparts in the Euclidean version of the model. 



A. Self-affine roughness 

Various statistical roughness estimators can be used 
to characterize the self-affine properties of a rough front. 
Consider for example the standard deviation a(Ax) of 
the height differences between points separated by a dis- 
tance Ax. A self-affine front obeys cr(Ax) oc Ax^ . Simi- 
larly the width w(L) of a front of length L (i.e. the stan- 
dard deviation of the height distribution along the front) 
scales as w(L) oc L**. Fourier or wavelet transforms are 
also of standard use. 
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FIG. 3: Width w(L) of the front (standard deviation of the 
height distribution) for growing lateral sizes L after normal- 
ization by a power law of exponent £. The simulations have 
been run from over 100 x 10 6 iterations steps for L = 2 s up to 
25 x 10 9 for L = 2 20 . The circle and the square symbols corre- 
spond respectively to a uniform and a Gaussian distribution 
of the trap depth. 

In the context of this study we can also design a 
"wavelet like" roughness estimator which exploits the 
natural tree structure associated to the ultrametric dis- 
tance. At the level £ = n = log 2 L we define co 2 (n) as the 
variance of the height difference between nearest neigh- 
bors: 

J{n) = Uh i -h j f) (6) 

\ / d u {i,jj— 2 n — 1 

At the upper level the height of a node is thus simply 
defined as the arithmetic average of its two ancestors: 
hl e -V(i) = [ftW(2i)-ftW(2i + l)]/2 and the correspond- 
ing variance w(m) is computed. This sequence is iterated 
up to the root of the tree. At each level m corresponds 
an ultrametric distance m = n — I + 1 and we have 

Lu 2 (m) cx 2 2C(m - 1} , or w(m) oc d u (m) c (7) 

We present now numerical results obtained for these 
various roughness estimators. The simulations have been 
performed on systems of sizes up to L = 2 20 . In Fig. |3| 
and 0| we show the scaling behaviors obtained for the 
wavelet roughness estimator u>(d u ) and the width w(L) 
of the interface. We obtain perfect power law behaviors 
over six decades and we only show here in logarithmic 
scale the residuals after normalization by a power law of 
exponent £. This procedure is a very sensitive way of de- 
tecting deviations from a power law. Note that previous 
published works deal with log 2 L < 10. 

We present simulations performed with two kinds of 
distributions for the trap depths, uniform and Gaussian 
(respectively denoted by the subscripts a and b in the 
figures). 

An estimation of the roughness exponent can be ex- 
tracted from each individual set of data. Note that 
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FIG. 4: Wavelet roughness estimator u)(m) against the ul- 
trametric distance d u (m) = 2 m — 1. The simulations have 
been performed on a system of size L = 2 20 run over 25 x 10 9 
iteration steps. The circle and the square symbols correspond 
respectively to a uniform and a Gaussian distribution of the 
trap depth. 



the fluctuations of this estimate due to the choice of 
the nature of disorder or the roughness estimator are 
larger than the deviation form the power law behavior it- 
self. All results are presented here with the central value 
C = 0.391. This value slightly underestimates the results 
obtained from the width estimator and slightly overesti- 
mates the results obtained from the wavelet estimator. 
A conservative estimate of the roughness exponent thus 
appears to be: 

C = 0.391 ± 0.005 (8) 



B. Avalanche behavior — Dynamic exponent 

The avalanche behavior is characteristic of the dynam- 
ics of the front close to threshold forcing. Although an 
extremal driving does not allow to recover the real dy- 
namics of the front, avalanches associated to a given level 
of the driving force F can be reconstructed from the his- 
tory of the extremal force signal s(t). An avalanche thus 
consists of a continuous series of depinning events such 
that s(t) < F. Instead of reconstructing these avalanches 
it is classical in the framework of extremal models [2{| 
to work directly from the extremal force signal. Let us 
consider a time lag At (a number of iterations), we intro- 
duce the distance Ax along the front between the sites 
depinning at t and t + At respectively. It appears @ that 
the distributions of these distances Ax at fixed time lad 
At can be rescaled on a universal form: 

P^AO-jL,*^) (9) 

where ipiu) oc u a for u <C 1 and ip(u) w cste for u> 1. 
The exponent a is well approximated by the exponent of 




FIG. 5: After rescaling the distribution of ultrametric dis- 
tances d u between sites depinning at time t and t+At collapse 
on a unique master curve. The dynamic exponent used for 
the rescaling is z — 1.39. 



FIG. 6: Distribution of the rescaled width of the depinning 
front for sizes L = 2 10 , 2 11 and L = 2 12 in the Euclidean 
model and L = 2 10 , 2 12 and L — 2 14 in the ultrametric 
model. All distributions collapse on the same curve. 



elastic kernel of Eq. (2) a f» 2 and the dynamic exponent 
z can be related to the roughness exponent: if a sequence 
of At depinning events spreads over a distance Ax along 
the front, the knowledge of the roughness of the front 
over Ax leads to At w AxAx c thus z = 1 + CUE!. 

This scaling is recovered in the framework of ultramet- 
ric depinning, where we measured P(d u , At). We check 
on Fig. [HJthat after rescaling all distributions collapse on 
a unique master curve. The rescaling was obtained with 
the value z = 1.39, for large arguments the behavior of tp 
is well approximated by a power law of exponent a = 2. 



V. SCALING FUNCTIONS 

In addition to critical exponents, the pinned state is 
also characterized by scaling functions 27] . In the specific 
case of interface dynamics, this means that the critical 
properties of the front are described by a universal expo- 
nent and a universal function describing the fluctuations 
of the (rescaled) width of the interface. This property 
has been evidenced for various growth models (Edwards- 
Wilkinson, Kardar-Parisi-Zhang. ..) [27] and has been re- 
cently applied to the case of depinning interfaces [28| . 
Note that in the latter case, beyond the interface width, 
the technique can be used to characterize other fluctuat- 
ing quantities. In particular, the distribution of the de- 
pinning threshold of a finite elastic line under extremal 
driving can be shown to be universal|23l |29j. In the 
present study we show that the choice of an ultramet- 
ric distance does not affect these universal distributions. 
More precisely the fluctuations of interface width and de- 
pinning threshold are shown to be described by universal 
functions and these functions appear to be very close to 
or identical to their counterparts obtained in the frame- 
work of a depinning model with Euclidean distance. 

These statistical distributions are however sensitive to 



the boundary conditions (periodic boundary conditions 
vs isolated system) 30]. In the following we use peri- 
odic boundary conditions. In the ultrametric case the re- 
summation over all replicas induces an additional mean 
field contribution 1/2L 2 equally shared by all Gij,i ^ j. 



A. Universal width fluctuations 

Following Ref. we study the distribution of the 
rescaled width w 2 / {w 2 } where (w 2 ) is the temporal aver- 
age of the width of a depinning front of finite extent L. 
In Fig. we present results obtained for various system 
sizes in both cases of Euclidean and ultrametric distance. 
We observe that all rescaled distributions collapse onto 
a master curve whatever the size of the front and the 
Euclidean or hierarchical metric. 



B. Universal depinning force fluctuations 

As developed above, the extremal dynamics gives a di- 
rect access to the fluctuations of the driving force needed 
to depin the front site by site. Most of these fluctua- 
tions simply correspond to the depinning of near neigh- 
bors (which receive the largest contributions of the elastic 
redistribution) and are highly sensitive to the details of 
the pinning force disorder. The other part of these fluc- 
tuations corresponds to depinning events taking place at 
a larger distance d from the previous depinning site. In 
other terms, the front can be regarded as pinned over 
the distance d between the two successive extremal sites. 
Conditioning the depinning force distribution to this dis- 
tance d between successive extremal sites thus allows to 
define distributions of size dependent effective thresholds. 
(Note that in the context of fracture, it can be seen as 
a distribution of effective toughness |3lj])- An effective 




FIG. 7: Distribution of depinning forces (bold line) and con- 
tributions corresponding to growing distances along the front 
between successive depinning sites. The larger the distance 
the narrower the distribution and the closer from the critical 
threshold s*. 

threshold over a distance d can also be seen as the force 
needed to entirely depin a front of size d. It can actually 
be shown [H H<| that as d increases, the distributions 
P(s; d) tend to peak and approach the critical threshold 
s* . This behavior is recovered in our ultrametric vari- 
ant as illustrated on Fig. [7| where we plotted the global 
distribution P(s) and the contributions P(s; d u ) corre- 
sponding to the different ultrametric distances. 

The typical elastic force fluctuations over a distance 
a can moreover be estimated: Sf el oc (H x ~ c ). Both 
the width <r s (d) and the gap to the critical threshold 
S s (d) = s* — (s)(d) of the center of these distribution 
actually follow this scaling. As developed in Ref. [2j| 
this property allows for a precise extrapolation of the crit- 
ical threshold by extrapolation of the linear relationship 
between (s)(d) and a s (d) up to the force value canceling 
a s . Moreover this rescaling results in a collapse of all 
distribution over a single master curve: 

P(s; d)=d 1 - i X [d 1 - C (s*-s)] . (10) 

We see on Fig. [3] that the choice of an ultrametric 
distance slightly changes the shape of the distribution 
obtained after rescaling. This may mean that these con- 
ditional distributions are more sensitive to the boundary 
conditions and the loss of translation invariance induced 
by the ultrametric model than a macroscopic quantity 
such as the width of the interface. 



VI. CONCLUSION 

We developed a depinning model based on the use of 
an ultrametric distance. This choice allows to reduce 
the complexity of an elementary step of computation to 
log 2 L with respect to the system size L. We performed 



d=2 3 (Euclidean) ° \ 




4 



(s-<s>)/a 



FIG. 8: After renormalization the centered unit conditional 
distribution of depinning force fall onto a single curve inde- 
pendently of the length d. The master curve obtained in the 
ultrametric case (continuous lines) slightly differ from its Eu- 
clidean counterpart (symbols). 

numerical simulations on systems of size 10 6 more than 
two orders of magnitude larger than in previously pub- 
lished results. We propose an estimate of the roughness 
exponent ( = 0.390 ± 0.01 consistent with the values 
measured for the Euclidean counterpart of the model. 
Moreover we obtained scaling functions either identical 
or very close in both cases. Therefore the choice of an ul- 
trametric distance appears not to affect the universality 
class of the depinning transition. 

Beyond the numerical efficiency, this model may thus 
serve as a starting point for a real space renormalization 
analysis. Indeed a coarse grained picture of the model 
preserves exactly the same structure with "dressed" 
thresholds and front advances, and thus the solution of 
a simple L = 2 front model may open the way to an an- 
alytic determination of critical exponents and universal 
scaling functions. 



APPENDIX A: NUMERICAL 
IMPLEMENTATION OF THE MODEL 

In this appendix we describe the algorithm used to run 
this ultrametric model of depinning. As above mentioned 
the structure of the algorithm is roughly similar to the 
original one but steps B and E are characterized by a 
complexity in log 2 L instead of L. 

Let us first consider the redistribution of elastic forces 
that take place after a depinning event. The increment 
of force at a given site j only depends on its distance 
d u (io,j) — 2 m — 1 to the extremal site io. Actually 2 m_1 
different sites are at this same distance from the extremal 
site io and form a subtree (see Fig. |2J). Instead of up- 
dating sequentially the elastic force increment on every 
sites it is possible to update the elastic force acting on 
the whole subtree. To this aim we define the force ipf t 
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as the force transported by the i branch of the £ level 
of the tree. The first level (£ = 1) here corresponds to 
the two branches at the root and the level (£ = log 2 L) 
to the branches pointing toward the L sites of the front 
(the leaves of the tree). Using this definition the elastic 
force acting on site i is nothing but the sum of the forces 
transported by its ancestor branches: 

n— log 2 L 

t=i 

In the hierarchical algorithm the step E thus con- 
sists of updating the branches of the force tree affected 
by the redistribution of elastic forces : ipj l Q t = ipj l Q e + 

Gi 0tm 5hi , jo = {*o/2 n } where {«} is the operation 
of exchanging the parity bit of i. 

In the same spirit we can build a tree to determine the 
value and location of the extremal site. The determina- 
tion of the extremal site requires a priori L — l operations 
of elementary comparisons to find s* — min^^ — /• ). 
After each depinning event, the elastic forces are updated 
on every sites so that the same L — l operations have to 
be performed at each iteration. In the ultrametric model 
the situation is slightly different because the elastic force 
on a site can be written as a sum of the forces trans- 
mitted by the branches of the tree (see Eq. IA1|) . As a 
consequence two sites i and j share all force components 
acting on branches above their nearest common ancestor. 



then we proceed iteratively up to the root of the tree: 



cr 4! p_i = min (a 2i ,p + <pf iiP , a 2 i+i, p + V>2i+l, P ) ( A6 ) 



a itP -i = a 2l , p if cr 2l ,p + tpf iiP < a 2l +i, P + f2i+i, P 

= OC2i+l,p if C2i,p + f2i,p > °2i+l,p + ffi+l.p 

(A7) 

The location and the value of the extremal site are thus 
given by 

s* = (71,0 , io = ai t o (A8) 

Without prior knowledge, the computation for the deter- 
mination of the extremal site of a given conformation is 
simply the sum of all comparisons at each level of the tree: 
Y^n 2^-1 _ 2" — 1 = L—l. This is exactly the same result 
as in the standard case. Consider now the situation after 
the depinning event: the trap depth f,o is updated at the 
extremal site and only n = log 2 L branches of the tree 
are altered by the elastic force distribution. A sequence 
of n comparisons thus allows to find the new extremal 
site: we start at the former extremal site io, 



p = n : j = io , a jo = t jo > <*j = *o (A9) 



ft = *, 



n— log 2 L 

n— log 2 L 
£=m+l 2 " 



± ■)■. 



(A2) 

(A3) 
(A4) 



We can then perform the comparison operations level 
by level going up the tree of elastic forces. To be more 
specific let us define the following hierarchical structure. 
At level n = log 2 L we define 



(A5) 



then for p = n to p = 1 we proceed iteratively to 
update the tree and determine the value c^o and location 
Qio of the new extremal site. 



cr 



j p =jo/2 n - p , k p = {j p } 



a j P -i,p~i — a j P ,p a jp-i,p-i ~ a j P ,p 

~ a j P ,P ^ a jp-i,P-l ~ a k p ,p 



(A10) 



(All) 



(A12) 
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